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Abstract 

This is the second paper on a new formalism for relativistic quantum measurements. 
Here, we construct a fully relativistic model for detectors that takes into account the de¬ 
tector’s state of motion, intrinsics dynamics, initial states and couplings to the measured 
field. The dual classical/quantum description of the detector is implemented by using a 
master-equation type of approximation for the coarse-grained pointer variables. Then 
we identify the probabilities that correspond to ideal measurements, i.e., measurements 
that are largely insensitive to modeling details of the apparatus. The Unruh-Dewitt 
and Glauber detectors are recovered at the appropriate limits. We employ our results 
to models of particle detection, photodetection and relativistic spin measurements, and 
we derive an ideal distribution for relativistic time-of-arrival probabilities. 


1 Introduction 

This is the second in a series of papers on relativistic quantum measurements, presenting 
a method that applies (i) to any Quantum Field Theory (QFT) and for any held-detector 
coupling, (ii) for the measurement of any observable, and (hi) for arbitrary size, shape and 
motion of the detector. 

The hrst paper of the series [1] contains the main ideas of the theory including the math¬ 
ematical set-up and the physical assumptions. Its main result was the derivation of the 
probabilities associated to any n distinct measurement events. These probabilities depend 
on a 2n-correlation function on the QFT, and also on the physical characteristics of the appa¬ 
ratuses, including their states of motion. All information about each apparatus is contained 
in a single function, the detector kernel. 

In this paper, we develop a fully relativistic model for detectors, and we derive a general 
expression for the detector kernel that is valid in any set-up. Then, we specialize to three 
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cases of interest: particle detection, photo-detection and relativistic spin measurements. We 
show that our theory recovers Glauber’s |2] and Unruh-Dewitt detectors [3li at appropriate 

limits. 

Throughout this work, we employ the word ’’detector” or ’’apparatus” to refer to a de¬ 
tecting element whose records can be correlated with a single microscopic event. We do not 
need to assume that the detecting element is elementary in the sense of [3], like for instance, a 
single bubble in a bubble chamber or a wire segment in a wire chamber. We assume that the 
detecting element records only a single event in each run of the particular experiment. In this 
sense, an experiment with n distinct records of observation, requires n distinct independent 
detecting elements, one for each record. 

The consideration of temporal coarse-graining is essential for the consistent dehnition of 
detection probabilities. In particular, temporal coarse-graining is necessary for a unique 
determination of the time associated to a measurement event. A record of observation is 
typically an excitation of hnite spatial extension. Different ’’parts” of the excitation move 
along different classical trajectories, so they correspond to different proper times. Unless one 
coarse-grains time at a scale a much larger than the size of the record, the attribution of a 
time instant to a measurement event is ambiguous. Thus, relativity sets a lower bound to 
the scale of temporal coarse-graining characterizing an apparatus. Furthermore, the coarse- 
graining time scale also depends on the intrinsic dynamics of the detector and on the effects 
of the environment. This dependence is incorporated into a function rj of the proper time r, 
which we call the degradation function. The degradation function describes the persistence 
of measurement records in time. 

Our modeling of the measurement apparatus leads to a detection kernel that is dehned 
through a specihcation of the following parameters. 

(i) an embedding function that describes a detector’s world tube in spacetime, i.e., it 
represents the detector’s macroscopic motion; 

(ii) a set of positive operators describing the measurement of variables other than the 
spacetime coordinates of the measurement event (e.g., spin, momentum); 

(iii) the spatial and temporal resolution of the spacetime coordinates; and 

(iv) the degradation function. 

Each choice of the above parameters corresponds to a physically distinct apparatus. In 
some regimes, the resulting probability distributions are largely insensitive to hne details in 
the modeling of the measurement apparatus. These regimes correspond to ideal measure¬ 
ments. We note that, with this dehnition, ideal measurements are not, in general, projective 
measurements. 

After deriving a general expression for the detection kernel, we apply the formalism to 
specihc types of measurement. First, we derive an ideal probability distribution for the time- 
of-arrival of relativistic particles. This probability distribution generalizes the relativistic 
POVM of Ref. [6] and is compatible with different couplings of the held to the measurement 
apparatus. We also consider the case of photodetection, where we show that Glauber’s theory 
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is obtained at the limit that a is much smaller than any time-scale characterizing the state of 
the electromagnetic (EM) held. Finally, we derive a general expression the time-integrated 
probability density associated to relativistic spin measurements. 

This paper is structured as follows. In Sec. 2, we briehy summarize the results of Ref. 
[T], setting the background for the present work. In Sec. 3, we present an important approx¬ 
imation for the self-dynamics of the internal degrees of freedom of the apparatus. In Sec. 4. 
we present our model for the measurement apparatus, and derive a workable expression for 
the detector kernel. This is employed in Sec. 5, in order to derive a general characterization 
of relativistic detectors. In the following three sections, we present applications of the for¬ 
malism: particle detection (Sec. 6), photo-detection (Sec. 7) and spin measurements (Sec. 
8). In Sec. 9, we summarize and discuss our results. 

2 Summary of past work 

The main result of Ref. pQ is the construction of a probability density associated to n 
measurement events. The physical system under consideration is described by a QFT on 
Minkowski spacetime. The quantum helds (priX) are dehned on a Hilbert space X and they 
transform covariantly under a representation of the Poincare group on T. The helds interact 
with n independent detectors, resulting to n independent measurement records. 

Spacetime description of detectors. Each detector admits a dual-quantum classical de¬ 
scription. The spatial extension of the detector at rest corresponds to a subset S of Euclidean 
space . We represent the points of S by three-vectors q. A moving detector corresponds 
to a world-tube S' x i? in Minkowski spacetime, described by a spacelike embedding function 


S :SxR^M. (1) 

The embedding is expressed in terms of coordinate functions £^{T,q), for r G /2 and 
q & S. Given the embedding function, we dehne the frame vector helds on the world-tube 

E^ = 

dr * dqi ■ 

The time-variable r is the proper time of the path E'^{T,q) for hxed q, i.e., so that 
T]fj_yE^E^ = — 1. We assume that each measurement event is specihed by an accuracy of order 
a in the proper time; a is the temporal coarse-graining scale for the apparatus. The, the 
detector’s size must be much smaller than a, so that a single proper time parameter suffices 
for the temporal ordering of recorded events. 

Quantum description of detectors. Each detector is associated to a Hilbert space /C that 
splits into two complementary subspaces: /C = KX © /C“. The subspace /C“ corresponds to 
the absence and the subspace to the presence of a macroscopic measurement record. We 
denote the projection operator associated to /C+ by E, and assume an initial state |a;) G /C_. 
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Consider, for example, an avalanche diode. K, represents the Hilbert space of holes and 
electrons in a jnnction, snbject to snitable bonndary conditions (external voltage, tempera- 
tnre and so on). Let I be an operator that corresponds to the cnrrent throngh the diode. 
Then, /C_ is the snbspace corresponding to all eigenvalnes of I smaller than some critical valne 
Iq (signihcantly larger than the cnrrent dne to thermal flnctnations) and /C+ corresponds to 
all eigenvalnes of I larger than Iq. 

The self-dynamics of a detector is described by a Hamiltonian operator h on /C. The 
Hamiltonian h generates translations with respect to the proper time r. The self-dynamics 
of the evolntion operator onght to leave the snbspaces /C± invariant, otherwise it wonld 
continnonsly generate measnrement records, even in absence of a measnred system. Hence, 
[E, h] = 0. 

The interaction Hamiltonian describing the conpling between held and detector is 


V = j d^xYA{x)J^{x) (3) 

where Ya{x) is a composite operator on IF that is a local fnnctional of the helds 0^, and 
J^{x) are cnrrent operators dehned on the Hilbert space K, of the detector. H is a collective 
index for the composite operators. 

Observables. Let A stand for all observables recorded by the detector. We denote the 
associated positive operators on /C by F{X) . In absence of a detection event, there are no 
observables so the operators F{X) are exhanstive on the snbspace associated to detection, 
Exnx)^E. 

The set of measnrement ontcomes always inclndes the location of the detection event. 
Hence, the parameter A in the positive operators F{X) is a shorthand for the pair 
where Q is a pointer variable that correlated to the coordinate q of the detector, and fi refers 
to pointer variables for magnitndes other than position (e.g., momentnm and spin). 

Both p and Q are highly coarse-grained variables, since they corresponds to macroscopic 
records. Hence, F{Q,fi) can be expressed as a prodnct 

F{Q,fr) = F,iQ)FM, (4) 

where Fi{Q) and T 2 (/i) are POVMs for the pointer variables Q and fi [7]. This implies that 
/ (PQF^iQ) = E and = E. 

Stationarity condition. We assnme that the initial state Iw) of the detector satishes the 
stationarity condition 

J^{x)e-^^^\u) = J d^qJ^{q)6%Y -S\q,T))\uj'), (5) 

for some vector |a;') and cnrrent operators J"^{q) dehned on S. The stationarity condition 
l|S|) implies that the cnrrent operators and the initial state are combined in snch a way snch 
that the only time dependence of J^(x, r)e“*^'^|(x;) is dne to the motion of the apparatns. 
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For a static detector, Eq. (Ej) means that |a;) is an eigenstate of h, and the corresponding 
energy eigenvalue is chosen to be zero. For a moving detector, Eq. ([5]) means that the 
Hamiltonian h affects only the part of the quantum state that corresponds to the apparatus’ 
macroscopic motion. 

Probabilities. Given the assumptions above, we compute the probability density for a 
sequence of n measurement events by n detectors. Each event is characterized by a value of 
the proper time variable r of the associated detector, one value for the observable Q that is 
correlated to the position, and other observables fi. It is important to emphasize that both 
the time and space coordinates are treated as random variables, and the probabilities dehne 
densities with respect to time. This is achieved using the Quantum Temporal Probabilities 
(QTP) method [B] that has been developed for addressing the time-of-arrival problem in 
quantum mechanics nn- 

In the QTP method, the time of an event is treated as a quasiclassical macroscopic vari¬ 
able. For such variables, probabilities can be rigorously dehned. This implies that we must 
introduce a coarse-graining time scale a that corresponds to the minimal temporal resolution 
of any event recorded by the detector. In the present context, the scale a refers to proper 
time rather than coordinate time. We shall explain the origins and physical interpretation 
of a in the next section. 

In Ref. [1], we computed the probability density associated to n detection events, each 
corresponding to a specihc value of the proper time r*, of the coordinate and the observable 
/ij of the Tth detector. 


Pipii Qli pit'^2-1 Q2-I P‘21 ■ ■ ■ I'T'niQni Pn) 

Y. f ■■■dXr, f d^X [... dX'^G{X,, Hi;..., X,,, Rp ...; X^, RJ 

A r> ID J 

xi?(i)(Xi, Hi; Xj, i?i|ri, Qi, /^i)... R(n){Xni Hi; X^, i?i|r„, Pn)i (6) 




where i = 1,2,... ,n labels the apparatuses. The correlation function G in Eq. ([B]) has 2n 
arguments and is dehned as 


Tr 


G(Xi, Hi;..., X„, H„|X;, Hi;...; X;, = 

riPjA'j.. .p,(A'i)|Ar|FB,(A")... Fb.(a;)| 


(7) 


where po is the initial state of the quantum held, X = {x,t) and Ya{X) = Ya{xR) = 
is the Heisenberg picture composite operator of the interaction Hamiltonian. 
We note that the correlation function G has n time-ordered arguments (T) and n arguments 
in reversed time order (T). 

All information about the Tthe detector is contained in the function Rp), the detector 
kernel, which is dehned as follows. 


R(X,A-,X',B\T,Q,fi) = j dsg.,(s) j tPq j Sq'S^X - Sir + ^-,q)]S‘^[X' - Sir 


s 

2 ’ 


a 

( 8 ) 
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In Eq. (IE]), f/tr is a function of proper time which can be conveniently chosen as a Gaussian 

g^{s) = exp[-sV(8a2)]- 

In the following two sections, we will elaborate on Eq. ([8]) and derive a simple equation 
for the detector kernel. 


3 The microscopic origins of the coarse-graining scale 


3.1 The master-equation approximation 


The probability density ([S]) is well dehned provided that the temporal variables are averaged 
at a time-scale a. The scale a corresponds to the irreducible temporal resolution of the 
apparatus. By ’’irreducible”, we mean that it arises out of the fundamental quantum de¬ 
scription of the apparatus. The actual resolution of an apparatus also includes contributions 
from classical uncertainties, for example, uncertainties in the propagation of the signal and 
uncertainties of the classical clock that determines the timing of an event. In what follows, 
we explain how a is inferred from the intrinsic dynamics of the apparatus. 

First, we recall von Neumann’s model of quantum measurements [S]. In this model, the 
self-Hamiltonian of the apparatus is negligible when compared to the interaction between the 
apparatus and the microscopic system. The interaction term causes the correlation between 
the microscopic observable and the pointer variable of the apparatus. One expects that a 
small but non-zero self-Hamiltonian induces fluctuations of the pointer variable around the 
recorded value, i.e., a degradation of the measurement signal. In our model, we provide a 
characterization of this degradation. 

Let us use the symbol A to denote all measurement records Q and fi in Eq. (|8]). The 
records are represented by the positive operators 11(A). They correspond to macroscopic 
pointer variables that are dehned in terms of collective degrees of freedom in the apparatus. 
This means that the operators F{X) describe highly coarse-grained observables. 

There are several different ways for quantifying the coarse-graining of a positive operator 
F(A). For example, we can dehne the associated density matrix 

F(A) 

^ TrF(A)’ 

and use the von Neumann entropy S\ = —Tr{px logp\) as a measure of coarse-graining. 

In Eq. (jH]), the detector’s self-Hamiltonian h appears in a term \/F(A)e“®^^\/F(A). We 
evaluate this product using the fact that F{X) is a highly coarse-grained operator. To this 
end, we adapt an approximation that is commonly employed in the derivation of master 
equations for macroscopic variables [12]. We write 

\/F(A)e-'^^\/F(A) = F{X)r]{X, s), (10) 


where the function 


r]{X,s) 


Tr 


F{X)e 


—ihs 


TrF(A) 


( 11 ) 
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will be referred to as the degradation function of the detectoi0. Eq. flTOl) is the master- 
equation approximation for the pointer variables. For analogous conditions of classicality for 
coarse-grained quantum variables, see Refs. 

For a justihcation of Eq. flTOl) . let us hrst consider the case of projective F{X), i.e., 

assume that ^ F{X) = F{X) and F{X)F{X') = 0, for A 7 ^ A'. For each A, F{X) projects into 
a subspace JCx of the Hilbert space /C+. 

Since A corresponds to macroscopically distinguishable observables, there is a concep¬ 
tual analogy between the subspaces }C\ and the macrostates of non-equilibrium statistical 
mechanics. In systems with a large number of degrees of freedom, the microscopic quan¬ 
tum dynamics are expected to generate an effective stochastic dynamics at the level of the 
macrostates. That is, we expect that the Hamiltonian evolution induces a probabilistic 
evolution between the subspaces K-x characterized by different values of A. The latter is ex¬ 
pressed in terms of a transition probability function G{X, 0; A', s) expressing the conditional 
probability that the recorded variable will take value A' at time s, provided it took value A 
initially. 

In non-equilibrium statistical mechanics, the transition probability G{X, 0; A', s) is com¬ 
monly estimated via a rerandomization condition, like Boltzmann’s Stosszahlansatz, or the 
transfer matrix approximation [13]. Here, we assume that for sufficiently large times s, the 
density matrix of the system is well approximated by the density matrix (|9|) of maximum 
ignorance on JCx- Then, 


^(A, 0; A', s) = Tr{F{X')e-^~^^pxFh 


Tr 


F(A')e-*'^*F(A)e 


■-.ihs 


TrF{X) 


( 12 ) 


Eq. flT^ is an analogue of Pauli’s master equation, in the sense that it monitors only the 
probabilities with respect to the different subspaces JCx and ignores cross-terms in the density 
matrix. 

The rerandomization assumption implies that the evolution operators are assumed to act 
only on states of the form px and produce states of the form px' with some multiplicative 
coefficient, i.e., it implies that 


F(A)e-“*F(A) = .,(A,s)F(A), 


(13) 


with p{X,s) given by Eq. (ITT]) . Substituting into Eq. flT^ . we hnd 

G(A,0;A,s) = |7^(A,s)|1 (14) 

Eq. (IT3|l was justihed for projective F{X). Continuity implies that its analogue also holds 
for approximate projectors, i.e., for F{X) that satisfy the conditions 


« 1 , 


Tr|F(A)F(A0| 

TrF(A) 


« 1 if A ^ A'. 


(15) 


Tr|F(A)-^F(A)| 

TrF(A) 

Then Eq. (ITnl) is obtained. 

^If the operators F{X) are not trace class, a suitable regularization is necessary before taking the quotient 
in Eq. (ITTl) . 
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3.2 Meaning of the degradation function 

It is easy to show that the degradation function satisfies 

l^(A,s)|<l (16) 

h(A,0) = l. (17) 

The degradation function r]{\, s) coincides with the persistence amplitude [lU UHl [16] 
of the density matrix px, Eq. ([9|). The square norm |77 (A,s)p is the associated persistence 
probability. We chose the characterization ’’degradation” rather than ’’persistence” for r 7 (A, s) 
because we want to emphasize the degradation of the record due to the intrinsic dynamics 
of the apparatus, rather than the persistence of some initial state. 

The squared norm of the degradation function is expected to be a decreasing function 
of s. This is a typical behavior for persistence probabilities |T6]; in fact, the decay is often 
exponential. This is consistent with the fact that the evolution generated by Eq. ca is 
typically diffusive, and ‘dilutes’ any initial concentration of the probability distribution at a 
specihc value of A. 

Let us denote by the time scale characterizing the decay of ri{X, s). Eq. (|H|) implies that 
if Td is much larger than the temporal coarse-graining time-scale a, the degradation function 
does not influence the probabilities: for s < a, ri{X,s) ~ v{X,0) = 1, while all contributions 
to the detector kernel from s » a vanish due to the function g„. In contrast, for a » Td, Pa 
is effectively unity. The key point is that the degradation function suppresses all interferences 
from amplitudes with detection times differing hj St » Td- Thus, choosing a of or larger 
guarantees the consistent definition of probabilities—see. Sec. 2.2 in Ref. [T]. 

We conclude that the optimum value for a that is compatible with the dehnition of 
probabilities and leads to maximal resolution in the time measurements is of the order of 
Td- Thus, the degradation function of the detector naturally determines the scale of temporal 
coarse-graining. 

For many observables, different values of the pointer variable do not correspond to differ¬ 
ent energies. For example, in a homogenous detector, there is no energy difference between 
detection records at different points. In such cases, the degradation function is insensitive to 
the specific value A, and we denote it as r]{s). 

In what follows, we estimate the degradation function and the associated timescale Td for 
two classes of measurements. 


3.3 Energy measurements. 


We consider an apparatus in which the recorded variable is the deposited energy, i.e., 
calorimeteiH. The positive operators associated to energy measurements are 


macroscopic record of energy absorption typically appears within an elementary detecting element, 
i.e., it corresponds to a transition in the energy levels of a single detecting element. Whether one focuses on 
the energy of a single element or of the whole detector is not an issue, because energy is an extensive variable. 
However, when trying to ascertain how the number N of particles in the detector is related to the decay scale 
Td, it is important to keep in mind that N refers to the number of particles of an elementary detector. 



( 18 ) 


#■(£)= fdE'f4E-E')pE,, 

where Pe are the spectral projectors of the Hamiltonian h. The sampling function is 

localized around E = Q with width A and satishes j dEfi\{E) = 1. 

Then, Eq. ffTTD gives 


, S dEU(E - E')e-‘’^'-g(E') 
JdEME-E'HE') 


(19) 


where g{E) is the density of states of the apparatus. 
We write g{E) = and we expand 


S{E') 


S{E) + 


1 

T(E) 


{E' - E) 


1 

2C{E)T{Ef 


{E' - Ef. 


( 20 ) 


The functions S{E), T ^{E) = dS/dE and C ^{E) = can be identihed with the entropy, 
the temperature and the heat capacity of the microcanonical distribution associated to h, 
respectively. 

For energies E » A, we can use a Gaussian sampling function, f/^{E) = 
to obtain 


ri{E, s) = exp 


2[A-2 + (GT2)-2] 


ii!|B+(T + ^ 


-1 


( 21 ) 


In the regime where the sampling accuracy is much larger than the thermal fluctuations of 
the energy, A >> y/CT, the degradation function simplihes 


ri{E, s) = e 2 


CT±^2 


( 22 ) 


The decay time-scale Td is 


VCT' 


(23) 


The heat capacity C typically scales with the number N of particles of the elementary detector 
and vanishes as T —)• 0. Hence, the decay time scale decreases with l/y/N and increases 
as T —)■ 0. 
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3.4 Position measurements. 

A typical record in a position measurement is the appearance of a localized excitation in an 
underlying medium, such as a fluid or a lattice. For example, wire chambers and semiconduc¬ 
tor detectors are characterized by electric current excitations, bubble chambers by excitations 
in collective degrees of freedom. The position observable Q corresponds to any variable that 
characterizes the location of the excitation. The size 6 of the excitation determines the 
accuracy of the position sampling. 

Degradation in a position measurement corresponds to the diffusion of the excitation in 
the underlying medium due to the random interaction between the medium’s components. 
In a homogeneous system, the degradation function rjiQ^s) should be the same for all points 
(5, so it does not depend on Q. 

At high temperatures, diffusion in a homogeneous system is well described by the standard 
Brownian motion. The excitation is described by a probability density p{q), where q is the 
position coordinate in the underlying medium, and the probability density evolves with the 
propagator 


(24) 

where D is the diffusion coefficient. 

A measurement record corresponds to a probability distribution psifl) located around 
9 = 0 with spread 5. In a homogeneous medium the location of the measurement event 
does not affect the degradation, so there is no loss of generality in choosing 9 = 0. The 
corresponding persistence probability is 

lh(s)P = / dWG'(9,0;9',s)p5(9)p5(g')- (25) 

The persistence probability fl2^ can be estimated by substituting (9 — 9 ')^ with the spread 
5^ in the propagator. Then, 




1 

(47rD|s|)^A 


4Z?|s| 


(26) 


Eq. fl26|) implies that the decay time Td is of the order of (5^/D, i.e., it is inversely 
proportional to the diffusion coefficient D. 


4 Evaluation of the detection kernel 

In this section, we continue with the calculation of the detector kernel and we identify a 
simple expression for its general form. We use this expression, in order to simplify the 
probabilities associated to n measurement records, and to identify regimes that correspond 
to ideal measurements. 
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4.1 Position localization 

Using the master equation approximation, the inner product 




( 27 ) 


that appears in the detector kernel, Eq. ([H]), becomes 




( 28 ) 


where we employed Eq. (jH). We assumed a homogeneous detector, so the degradation 
function for a position measurement is Q-independent. 

The currents define the interaction between the detector and the quantum field. The 

vector j^(g)|a;') is typically characterized by a microscopic length scale I that corresponds 
to the range of the interaction. For example, ^ may be the atomic scale (Bohr’s radius) or 
the nuclear scale. Whatever the value of £, it is much smaller than the macroscopic scale 6 
that characterizes the localization of a position record. This means that 

- for \Q -q\ < 5, Fi{Q)J^{q)\uj') ~ J^{q)\uj'), and 

- for IQ -q\ » S, Fi{Q)J^{q)\u') ~ 0. 

Thus, we can write 

F,iQ)J{q)\oj) = usiQ - q)Jiq)\u;), (29) 

in terms of a sampling function us : S ^ that satisfies J d?Qus{Q) = 1. 

Thus, the inner product (l28|) involves a product \/us{Q — q)us{Q — q'). We choose a 
Gaussian sampling function 

“ (27r^2)3/2^ ^ 

Then, 

\/UsiQ - q)u5{Q - q') = us(^- W5{q-q'). ( 31 ) 

where 


ws{q) = ( 32 ) 

Assuming a scale of observation for Q much larger than 5 we can approximate us{Q — q) 
with a delta function 6^{Q — q) —here q = {q + q')/2 . Writing r = q — q', the inner product 
fl2Sll becomes 

W = ,(f,, i,)i3[Q - (4 + 4')/2|«>i(4 - 4')(V| j'’(4')f'2(A) jy?)!!*)'). (33) 
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4.2 Matrix elements of observables 

To further simplify the inner product fl33|l . we consider the Hilbert space V of test functions 
for the current operators V is the space of square integrable functions Ja ■ S ^ C, 

with an inner product 

(fig) = / d-hfl{q)9A{q)- (34) 

We employ a Dirac-like notation and we write the vectors of V as j-). The elements \q,A) 
dehne generalized eigenstates of the position operator on S. 

Then, we express the inner product {q')F 2 {iJi)J^{q)\oj') on /C as an inner product 

on V, 

= C(q', B\f(iJ.)\q,A), (35) 

where C is a positive constant and T(/i) are positive operators on V. The constant C is 
determined as follows. 

We assume that all matrix elements of the operator -J^q) between states of /C_ vanish, 
because otherwise there would be no meaningful correlations between microscopic particles 
and detection records. Hence, 


Since = E, 


(i - E)J-\q)\Lj) = 0 


(36) 


Cj2{<l',B\f(A\q,A) = {u,'\P(q')j*{q)\A}. (37) 

The correlation function is expected to vanish for \q—q'\ » i, where i is the characteristic 
scale of the interaction. Hence, at scale much larger than £, it is proportional to a delta 
function 


{u\J^iq')J^{q)\u) ^ C5%q-q'h^^ = Ciq',B\q,A), (38) 

where 7 "^^ ( 9 ) is a tensor that plays the role of the metric for the discrete degrees of freedom. 
In most cases 7 ^'® will be the Kronecker delta 5^^. However, when the collective index A 
includes spacetime or spinor indices, 7 ^^ may be a non-trivial function of q. 

Eq. fl35]) serves as a dehnition for C. It also implies that the positive operators T(/r) 
dehne a POVM on V, since by Eq. fl57|) . 

Then, the inner product (jSS]) becomes 

W = Crjif,, s)S^lQ -(q + q')/2]ws(q - q')(q', Bjf(/ 2 )jq, A). (39) 


12 


4.3 General form of the detector kernel 

From Eqs. ([8]) and fl3^ . we obtain the following expression for the apparatns kernel, 

R{X, A]X\B\t,Q, ^) = C J dsga-{s)T]{n,s) J d^rws{r) 

-S(r+^-,Q + - £(r - B\f{^,)\Q + 4), (40) 

where we set r — q — q'. 

When snmming over all valnes of the observed variable /r, a term of the form 

J2W,B\np)\i!,A)v(.P,s) (41) 

appears in the kernel. This term simplihes signihcantly when the degradation fnnction s) 
is /x-independent. Then, the snm of Eq. fjdTll eqnals 5^{q' — q)'j^^iq)r]{s), and we obtain 

R(X, /I; A", B\t,Q,ii) = j dsC-f'”g,(s)r,(s)S*{X - f (r + Q)]i‘[A' - f (t - |,Q)](42) 

The parameter C involves the conpling constants of the interaction between held and 
detectors and it is proportional to the total nnmber of particles that left a detection record. 
When we normalize the probabilities by dividing with the fraction of the particles that have 
been detected C drops ont. 

The detector kernel has a strong dependence on the temporal coarse-graining scale a and 
on the spatial coarse-graining scale d. A measnring event corresponds to a fonr-volnme of 
order in E x i2. The corresponding spacetime fonr-volnme depends on the detector’s 
motion as encoded in the embedding fnnction. We also recall that a >> 5, so that a single 
proper time parameter can be nsed for labeling the events in one detector. 

To snmmarize, we fonnd that each measnring apparatns is characterized by 

(i) its macroscopic motion in space-time encoded in the embedding fnnction T, 

(ii) a POVM T(/x) that is dehned on the Hilbert space V, describing the measnrement of 
all observables other than position, 

(iii) spatial resolntion d and temporal resolntion a of the spacetime coordinates, and 

(iv) the degradation fnnction ri{fi,T). 

5 Detection probabilities 

Next, we insert Eq. fl40|) for the apparatns kernel into the probability densities (jb]), and 
identify idealized expression for the probabilities associated to n measurement events of 
general validity. We hrst consider n = 1 and then proceed to the most general case. 


13 


5.1 A single measurement event 

We construct the probability density that the detector records an event charac¬ 

terized by proper time r, spatial coordinate Q and observable /r, 


= C j d^rdsg„{s)T]{ii,s)w5{r)GAB[S{T+ -^)] 

x(Q-^,S|f(/i)|Q + ^,A), (43) 

By Gab{.X,X') we denote the 2-point function ([7]) 


Gab(X,X') 


G(X,AlX',B)=Tr Va(X)poVb(X') 


(44) 


For observables p such that s) is p-independent, the integrated probability density is 
P{'r,Q) ■='^P{r,Q,p) = G I dsga{s)ri{s)G[S{T + ^,Q)]S{t - ^,Q)], (45) 

where G{X,X') = 7"^^Gab(W, X'). The same expression applies for the probability density 
for measurements with no records other than position, i.e, for F 2 {p) = E. 

The notion of an ideal measurement is of particular signihcance in quantum measurement 
theory. Ideal measurements are described by probabilities that are insensitive to features of 
the measurement device. For example, the most general measurement of position is described 
in terms of a POVM that commutes with the position operator. In an ideal measurement, 
the POVM consists of the spectral projectors of the position operator. In this case, it does 
not depend on specihc properties of the measuring apparatus, as for example, the localization 
scale of the measurement event. 

The probabilities depend on several different parameters with different physical sig¬ 
nihcance. Hence, there are several different types of ideal measurement. 


1. The point-like detector. In this regime, the detector’s world-tube collapses into a single 
time-like path, dehned, say, by Q = 0. Writing Xq{t) = T^(r, 0), we obtain 

F(r, p) = G j dsg^{s)ri{p, s)Gab[Xo{t + |), Xo(r - (46) 

where the matrices are dehned as 

= (47) 

In a point-like detector only ultra-local observables can be measured, i.e., observables 
that do not depend on derivative operators on the Hilbert space V. An example of 
point-like detector model is the Unruh-Dewitt detector Ell]. 
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2. The fully coherent detector. If the characteristic scales of the functions g^, ws and rj 

are much larger than any microscopic scales appearing in the correlation function Gab, 
but much smaller than the macroscopic scale of observation, we can set ws ^ l,ga ^ ^ 
and ri{fi, s) = 0) = 1, to obtain 

P{t,Q,i,) = C j ds j <iVG|£(T + |,0 + |)i£(r-|.e-^)] 

x((3-+(48) 

3. The incoherent detector. This condition applies only to apparatuses that record no 
observable /i, so that Eq. fH5|) applies. In this regime, the temporal coarse-graining 
parameter a is much smaller than any temporal parameter that characterizes the cor¬ 
relation function Gab- This is possible when the characteristic decay time-scale in r]{s) 
is itself much smaller than any temporal parameter in Gab- 

Then, we substitute g^ by a delta function ga{s) —?• '/^a6{s). Hence, 


P{r,Q,A = C'Gl£{Q,T),£{Q,T)] (49) 

where G' = y/^aG. The detection probability depends only on the diagonal elements 
of the correlation function. In Sec. 7, we show that Glauber photo-detectors [2] are 
special cases of incoherent detectors. 


5.2 Time-integrated probabilities 

When the time instant of a measurement event is not recorded, we employ the time-integrated 
probability densities 

roo 

P{Q,IX)= / dTP{T,Q,ft). (50) 

Jo 

The range of time-integration in Eq. fl50l) is [0, cxd), since we assume that the initial prepa¬ 
ration of the system corresponds to r = 0. Nonetheless, the probability density P(r, Q, g) 
can be analytically extended to negative values of r. In some set-ups, P{T,Q,g) is strongly 
suppressed for r < 0. This is the case, for example, in time-of-arrival measurements: par¬ 
ticles are prepared in a spatially localized state with support on momenta directed towards 
the detector. In such cases, we can extend the range of integration in Eq. fl50|) to (—oo, cxd). 

In two regimes, the time-integrated probability Eq. fl5U]l takes a particularly simple form. 
First, for the coherent detector, Eq. and for r integrated over (—cx3, cx)). 


P{Q,g) = G 


d^rTr 


Va(Q + \)mVB(Q - 




( 51 ) 
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where 


/ OO 

dTYA[£{T,Q)], 

■OO 


( 52 ) 


is the time-average of the composite operator along the detector’s world volume. 

Eq. flSTj) is a consequence of Eq. ([6]). However, the latter equation follows from an 
approximation by which we ignore a further classical smearing via convolution of the prob¬ 
ability densities with respect to time, at a scale of a [T]. Ignoring this smearing may lead 
to problems in the dehnability if the integral over r in Eq. fl5^ extends over the whole real 
axis. Without this approximation, the time-average composite operator should be expressed 
as 

/ CXD _ 

<iryU?)Pl£(T,(3)], (53) 

•OO 


where /o- is a smearing function at the scale of a. 

The second regime where Eq. (ISO]) corresponds to the incoherent detector. By Eq. fl49p . 


F(Q,/x) = C'7 


AB 


drTr 


Yb(S{t,Q))p„Yb(S{t,Q)) 


(54) 


Eqs. fl^ and (l5T|) apply to different physical situations. In Eq. (15T|) all superpositions 
associated to different instants of detection are decohered, while in Eq. flSTj) they are pre¬ 
served. Thus, the two equations may lead to different predictions in physical conhgurations 
that are characterized by very small temporal parameters. This is the case for for particle 
oscillation experiments [HI |T9] in which time-evolution induces a very small time due to 
the interference of particle states with different mass. In Ref. |6], it was shown that the 
standard formula for particle oscillations requires the postulate of an incoherent detector as 
in Eq. flM|) . In contrast, the consideration of a coherent detector as in Eq. flHT]) leads to a 
non-standard oscillation formula |2Uj . 


5.3 Multiple measurement events 

Next, we evaluate the probability density (E]) associated to n detection events, each corre¬ 
sponding to a specihc value of the proper time r*, of the coordinate and the observable 
/ij of the i-th detector. Each detector is characterized by an instantaneous space Si, and is 
spacetime motion is given by embedding function Si : Si x R ^ M. The temporal coarse- 
graining scale of the Tth detector is denoted as cXi, the spatial coarse-graining scale as 6i. 
and the associated the degradation function as rji^Ti). 
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We find 


= / Y[{CidSid^ri) 

^ i=l 
n 

i=l 

xG\£i{ti + —yQi + —), Ai; ...; f„(r„ + — 

fi(n - l-.O. - j), B.;...; £„(r„ -'^,Qn + jhB„], (55) 

where G(Wi, Ai;...; 74„|X{, Pi;...; X'^, B^) is the 2n-point function ([7]). 

In absence of records /ij, or for records /i* such that rji is /ij-independent, we obtain the 
integrated probabilities 


/ n n 

n n 

u.l,...,Un *=1 * = 1 


gaAsi)vi{si)] 


G|Bi(ti + j.Oi), -4i;... ;f„(T„ + - y,gi),^i;...; £„(r„ - B„],(56) 


As in the case of the single measurement, we can identify different regimes for ideal mea¬ 
surements, namely, point-like, coherent and incoherent detectors. They follow from Eq. fl5^ 
by an obvious application of the conditions presented in Sec. 3.3.1; their explicit expression 
is not given here. 


6 Particle detection 

We consider the probability density P{t,Q), Eq. fH5|) for the detection of relativistic particles 
of mass m, ignoring all other measurement records /r. Thus, the only observables are the 
spacetime coordinates for a particle-detection event. This set-up generalizes the treatment 
of time-of-arrival measurements of Ref. [6] by (i) considering different types of interaction 
and (ii) taking into account the detector’s spacetime trajectories. 

For simplicity, we assume spinless particles. The Hilbert space associated to a single 
particle is "Hi = L'^[R^, d^{k)] where djjt,{k) = 2cfn)^ujk Lorentz invariant integration 

measure. We consider a free scalar held 4>{X), so the held Hilbert space X is the bosonic 
Fock space associated to Pi. The held (f>{X) is dehned terms of the standard creation and 
annihilation operator on X: 

(j){X) = J d^i{k) , (57) 

where k ■ X = k^ = {uk,k), ojk = The creation and annihilation operators 

satisfy the commutation relations 

[dk, dk'] = [d\, d\,] = 0 [ofe, dl,] = {27rf{2ujk)6^{k - k'). (58) 
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6.1 Linear coupling 

First, we consider the case of a linear field-apparatns conpling, that is, we identify the 
composite operator Y{X) with (j){X). Since the interaction is linear to the creation and 
annihilation operators, this conpling corresponds to detection by particle absorption. 

We obtain 


= j dsr]{s){^\^[£{T - ^,Q)]^[£{t + (59) 

where we dropped the mnltiplicative constant C and absorbed into r]{s). In what follows, 
we assnme that ri{s) = r]{—s). 

We consider an initial state with a dehnite nnmber N of particles. Then, Eq. fl59l) 
becomes 


Pir,Q) = j dsr]{s)X^[8{T + |,( 5 ), - ^,Q)] 

+2N j dsriis) j (60) 

where 

A+(X,X') = (O|0(X')0(X)|O) = y (61) 

is the held’s Wightman’s fnnction and 

Piik, k') = y(T|aiafc/|T), (62) 

is the one-particle rednced density matrix associated to the state IT). The rednced density 
matrix pi is normalized as / dp{k)pi{k, k) = 1. 

The hrst term in Eq. fl60l) is state-independent and describes particle creation from 
the vacnnm as in the Unrnh effect. It has been extensively stndied in Ref. im. This 
terms vanishes for motion along inertial trajectories, and it is snppressed exponentially with 
increasing mass m. In general, its contribntion is very small, nnless the detector’s motion 
is characterized by time-scales mnch smaller than the coarse-graining time scale a. We will 
ignore this term in what follows. 

The second term in Eq. (I60D depends only on the one-particle rednced density matrix. 
This means that no correlations between particles can be accessed from a single local mea- 
snrement. 

6.2 The ideal time-of-arrival distribution 

We consider a static detector located far away from a particle sonrce, so that only particles 
with momentnm along the sonrce-detector axis can be recorded. Then, we can restrict to a 
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single spatial dimension. The corresponding embedding function is {£^,8^) = {t,Q), where 
Q > 0 denotes the distance of the detector from the particle source. Then, Eq. fl60|l becomes 

F(r, Q) = 2N J d^l{k)d^l{k')pi{k', k)p 


where r] is the Fourier transform of r]{s). Since we work in one spatial dimension, dfi{k) = 

dk 

(27r)(2(^fc) ■ 

We consider the time-integrated probability density P{Q) = drP^T, Q). For an initial 
state with support only on positive values of fc, and localized a.t x < Q, we can extend the 
range of integration for r to (— 00 , 00 ). Then 

P(Q) = N r dMk)pi(k, 4)^ (64) 

Jo \k\ 

Eq. fIMl) implies that fi{uk)/\k\ is proportional to the absorption coefficient a{u) in the 
rest frame of the detector. The absorption coefficient is defined as the fraction of particles 
with energy u absorbed per unit length at the detector. Hence, 

i)(a;) = Aa{oj)\/oj'^ — m^, (65) 


for some constant A. Eq. fl65l) imply that the degradation function for a detector can be 
reconstructed from the measurement of its absorption coefficient, i.e., from information that 
is macroscopically accessible. 

An idealized limit for this measurement corresponds to a constant absorption coefficient. 
We also consider a constant value of Q and normalize the probability distribution P(r, Q) as 
kT^d'rP{r,Q) = 1. Then, 


F(r,Q) = J dp{k)dp{k')pi{k, k')^J{uk + oJk'Y - ^gg^ 

We consider an initial pure state pi{k, k') = 'ipo{k)'ipQ{k') that is well localized in momen¬ 
tum, so that we can approximate {uk + oJk'Y — 4m^ ~ 4|fc||/c'|. Then, Eq. fl66|) becomes 


P{pQ) 


oJk 


(67) 


Eq. (|6711 is the POVM derived in Ref. |6]. It generalizes Kijowski’s POVM for the time 
of arrival [21] to relativistic particles. 


6.3 Moving detector 

For a general detector trajectory, we expand the embedding function r + s) around r 

and keep terms to first order in s, 

8^{Q,t + s) = 8 ^{Q,t) + u^{t)s, ( 68 ) 
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where = £^. Eq. fl68l) is a good approximation if aa « 1, where a is the proper 
acceleration of the path. 

Then, Eq. fl60|) becomes 


P{t,Q) j dfi{k)dfi{k')pi{k',k)fj ^i{k'-ky£ir,Q)^ 

i.e., fj is evaluated for energies in the rest frame of the detector. In the regime of the ideal 
detector of Eq. (l66ll . Eq. (j6^ becomes 

P{t,Q) = j dp,{k)dp,{k')pi{k', k)\/[{k + k') ■ u{t)Y — (70) 

Eq. (1701) dehnes the ideal probability distribution for the time-of-arrival in a moving detector. 


6.4 Quadratic coupling 


Next, we consider the case of quadratic held-apparatus coupling, that is, we identify the 
composite operator E(X) with the normal-ordered held product : :. The relevant 

interaction term involves one creation and one annihilation operator of the held and it cor¬ 
responds to detection through particle scattering. 

We consider a single-particle initial state for the held 

I'l'} = J dp(k)y!jo(k)dllO}. (71) 

Then, Eq. fjTHl) yields 

P(t,Q) = 2 j dsp(s) (^A+[£(r + ^,Q), (S(t - |,<3)]) 

+8 j dsp(s) j 


The hrst term in Eq. fl72l) is state-independent and describes particle creation from the 
vacuum as in the Unruh ehect. This terms is suppressed exponentially with increasing mass 
m. In general, its contribution is very small, unless the detector’s motion is characterized by 
time-scales of the order of the coarse-graining time scale a or smaller. We shall ignore it in 
what follows. 

For a detector on an inertial motion, the Wightman function Eq. fl72|) becomes 


A+(s) 


47r2 Jq Uk 


(73) 


We dehne k(s) := 477 ( 5 )A’''(s). Then, we restrict to one spatial dimension, with an 
embedding function (£^,£^) = (t,Q). Then, Eq. (I60|) becomes 


^k + ^k' 


i{uJk-ujy)T+i{k'-k)Q 
^ 1 


P('TiQ) = ‘^ / dp(k)dp(k')pi(k',k)K 


2 


(74) 







where k is the Fourier transform of k. 

Eq. frM|) coincides with Eq. fl63|) modulo the identihcation of k with Nfj. This implies 
that the quadratic held-detector coupling also leads to the ideal distribution fl66|l . 


7 Photo-detection 


Next, we apply our formalism to the case of photo-detection. The measured system is the 
quantum electromagnetic held 


Af,{X) = / dfi{k)'^\dk,r(^^,{k,r)e 


r=l 


Ak-X 




(fc,r)e 


-ik-X 


(75) 


expressed in terms of the standard creation and annihilation operators; e^{k, 1) and e^{k,2) 
are two orthogonal and transverse real-valued polarization vectors. 


7.1 Detection probabilities 

We consider an interaction between the electromagnetic held and the measuring apparatus 
of the form 


V 


SxF^yix) (g) 


(76) 


where = d^Ay—dpA^ is the held-strength tensor, and is the relativistic magnetization- 
polarization tensor associated to the measuring apparatus. 

The electric polarization is typically much larger than the magnetization in the rest frame 
of the detector. Thus, the dominant contribution to the held-detector coupling corresponds 
to a current p*(g) that represents polarization density in the rest frame of the detector. In 
this case, the stationarity condition, Eq. ([5]) takes the form 

M'^"(x)e-'^"|ca) = j d^q p'(g)<5^(a;' - T'(g, f))|a;'), (77) 


for some state vector |a;'). 

Following the reasoning of Secs. 3 and 4, we derive the (unnormalized) probability density 
P{t,Q) of detection at time r and position Q, 


PA,Q) = j dsgAs)v{s)G[£{Q,T + - |)]. 

The correlation function 


G{X, X') 




E,X')poE,{X) 




(78) 


(79) 
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is defined in terms of the electric field vector in the detector’s rest frame. 

Ei = \ {£“£•; - £‘'£f) F,. ( 80 ) 

the inverse of the indnced metric hij = on the instantaneons space S', and the 

initial state of the held po- 

Eq. (1781) applies to any photo-detector, irrespective of its state of motion. It is important 
to note that the electric held depends on Q and r not only throngh its argnment, bnt also 
throngh the tensor | j that projects onto the instantaneons space. 


7.2 Static detectors 


Next, we specialize to static detectors, characterized by an embedding fnnction S{t,Q) = 
{t,Q). In this case, the tensor does not depend on Q and r, and hij = 6ij . The 

electric held becomes 


/ z 

dn{k)uk ^ - al j{k, r)e 

r=l 


-ik-X 


( 81 ) 


where ujk = £^k^ = 1^1 and 

fi{k,T) = {S^S’^-rsnkMk,r), (82) 

are the nnit 3-vectors that describe polarization. 

It is convenient to proceed by introdncing the coherent states \z) of the electromagnetic 
held. They are dehned as eigenvalnes of the annihilation operator 


hk,r\^) = ^r{k)\z), (83) 

in terms of a sqnare-integrable complex-valned fnnction Zr{k). 

We employ the P representation for the initial state 

P = J Dzfp{z)\z){z\, (84) 

where fp{z) is a fnnctional known as the P-symbol of qnantnm state p [22]. The fnnctional 
integration in Eq. flMD is dehned in terms of the Ganssian measnre that dehnes the Bargmann 
representation of the qnantnm held [23] . 

The probability density Eq. (ITHD is expressed as 

P(r,Q) = j Dzfp{z)P,iT,Q), (85) 

where Pz{t,Q is the probability density Eq. fl78|) evalnated for a coherent state. Snbstitnting 
in Eq. flTHD . we hnd that P^ is a snm of three contribntions 

PziT,Q) = P^°\t,Q) + P^^^\t,Q) + P^^\t,Q), (86) 


22 



where 


P<°Ut,Q) 

P'.'\r,Q) 

P^'Ht.Q) 


j dix{k)r]{ujk) 

2 J dn{k)dn{k')uJkOJk'fj 

—2 / d^{k)dfi{k')ukUik'V 


+ ^k' 

2 

^k — ^k' 


Ck' ■ CfcC 


i{k~k')-Q—i{Ljf^—LUy)7 


Re 


Cfc ■ Cfc'C 




■ ( 87 ) 


We wrote Ck = X]r=i ^r{k)f{k, r), and we absorbed the term into the degradation fnnction 
g{s). 

The meaning of the three terms in Eq. fISB]) is the following. 


• The term is state-independent. It corresponds to a backgronnd noise of ‘false 
alarms’ in the detector [5] and its contribntion to the detection probability is negligible. 

• The term is generated by the terms in the density matrix that corresponds to a 
hxed nnmber of photons, i.e., by the restriction of the density matrix in the N photon 
snbspaces. It has the same form with the time-of-arrival probability distribntion fl55D 
for m = 0. 

• The term Pi corresponds to snperpositions in the nnmber of photons. It will be 
denoted as P 2 (pQ). 


7.3 The Rotating Wave Approximation and the Glauber detector 

The second term in Eq. 0861) involves a donble momentnm integral over (co¬ 
rotating terms), while the third term involves a donble integral over (connter- 

rotating terms). The dominant contribntions to the donble integral arise from valnes of k and 
k' that correspond to a slowly varying phase. This implies that co-rotating terms typically 
dominate over connter-rotating ones. It is therefore meaningfnl to invoke the Rotating Wave 
Approximation (RWA) and drop the contribntion of the connter rotating terms. 

The RWA is a dehning approximation of Glanber’s theory of photo-detection. The domain 
of validity of the RWA in photodetection has been a matter of some debate. Some anthors 
have shown that photodetection probabilities that are obtained with the RWA appear to 
violate cansality at short times and they proposed corrections to Glanber’s theory [2l]. Others 
indicate that a different form of the RWA in photodetection theory can gnarantee cansality 
[25] . A fundamental problem of the RWA is that it corresponds to a Hamiltonian that has 
no lower bonnd, and thns it is thermodynamically nnacceptable [2B]. However, it has been 
snggested that the proper renormalization of the interaction between detector and qnantnm 
held may lead to an effective RWA dynamics [27]. Regarding the degree of accnracy of the 
RWA and related approximations, see. Ref. [28| and references therein. 
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To examine the validity of the RWA, we consider a coherent state that corresponds to a 
pulse of mean wave-number ko localized around x = 0 at f = 0. We choose a Gaussian 

Cfc = (27r)3Co(27rA2)-3/2e-^, (88) 

where Co is a constant complex valued vector and A is the spectral width of the pulse. 

First, we evaluate the term Pz^\t,Q) in Eq. (1861) . We employ a saddle point approxi¬ 
mation. We expand ujk = ^ko +'Vko ' ~ ^q), where Vk = k/\k\ is the three- velocity vector 

for a photon. Then, we assume that fj varies slowly with u, so that the contribution of the 
degradation function is a constant eta{uko), we hnd 

(89) 

where Vk = k/\k\ is the three- velocity vector for a photon. 

We also evaluate the term P2{t,Q) in the saddle-point approximation. We hnd 

1 (k — k \ 

P^^\t,Q) = --e-^^Re Ij dsr]{s)e-^"^" j . (90) 

The key point is that the oscillatory terms in Eq. fl90l) must be averaged over a spatial 
region of size 6^ and for times of order a » 6. For Gaussian smearing functions, this 
averaging leads to suppression factors of order Thus, the contribution of the 

term Pz to the probability density fl86|) is suppressed, and the RWA is justihed. We expect 
that the suppression of the counter-rotating terms due to spatial and temporal coarse-graining 
in the detector holds for a large class of initial states. 

The RWA is equivalent with the rewriting of Eq. (1781) as 

= j dsri{s)GRWA[S{r + - ^,Q)]. ( 91 ) 

Eq. fl?IH involves a different correlation function 

GnwAiX.X') = 2W(vI/|4(-)(A')eW(A)|vI/), (92) 

where and E^~'^ denote the positive and negative frequency part of Ei respectively. 

For an incoherent detector, the probability distribution P{t,Q) becomes 

P{pQ) — Grwa[£-{TiQ)i£-{TiQ)]i (93) 

and, thus, coincides with the standard formula of Glauber’s photo-detection theory. 

For a single measurement event we recover Glauber’s photo-detection theory by assuming 
(i) the RWA and (ii) the incoherent detector limit. As we showed, the RWA is not an 
independent assumption, as it arises due to the spatial and temporal coarse-graining of 
the detector for a large class of initial states. Thus, the defining assumption of Glauber’s 
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theory is that the temporal coarse-graining scale is much smaller than any time parameter 
characterizing the initial state of the held. 

In Glauber’s theory, the n-time probabilities 

F(ri, Qi, /ii;...; r„, = Grwa[^ (n, Qi),..., S(Tr„QJ; ..., £ (r„, Q^)] (94) 

are dehned in terms of correlation functions 

inh .. .i-J-Tr (B<y(X„)... • (95) 

Eq. fl94)) follows from Eq. fl5^ assuming the RWA and the incoherent detector limits. 
However, we are skeptical about the validity of the RWA for general set-ups of n measure¬ 
ments. The RWA may misrepresent the effects of retarded propagation of the electromagnetic 
held |2n], hence, it may fail in conhgurations that involve several several detection events that 
are spacelike separated. In contrast, Eq. flSB]l is applicable in such systems. 


8 Spin measurements 


Our last example involves the calculation of probabilities for spin measurements in relativis¬ 
tic electrons. This issue is of particular interest for relativistic quantum information because 
there are diherent proposals about the operator representation of spin measurements in rel¬ 
ativistic systems [301 El] • 

Spin measurements are often described in terms of the Pauli-Lubanski vector |32] and 
related operator. This description is covariant, in the sense that the expectation values for 
the spin measurements are the same in all reference frames. However, as noted in Ref. |33] . 
there is no known physical coupling of relativistic particles to the detector compatible with 
this description. In fact, a careful description of spin measurements (as, for example, in 
Stern-Gerlach-type experiments [M]) is important in order to avoid violations of relativistic 
causality [35] . 

Our starting point in this section is the QFT description of particles of mass m and spin 
s = |. We consider the free Dirac held. 


^{X) = 'ip^+\X)+^^-\X), i{X) = i^+\X) + i)^-\X) (96) 
2 „ ^ 

^/,(+)(x) = ^ dij{k)cr{k)u{k,r)e^^-^ = J2 dn{k)dl{k)v{k,r)e-^'^-^ (97) 


r=l 

2 


r=l 

2 


'ijj^~^\X) =/ dfi{k)dr{k)v{k,r)e’''^'^ V’*' / dfi{k)cl{k)u{k,r)e (98) 


r=l 


r=l 


where Cr{k) and dr{k) are annihilation operators for electrons and positrons respectively. The 
four-vectors k satisfy = —m?. The Dirac spinors u{k,r) and v{k,r) are normalized as 

u{k,r)u{k,r) = 2m6rr' u\k,r)u{k,r) = 2ukSrr', (99) 

v{k,r)v{k,r) = —2m6rr' {k,r)v{k,r) = 2ukSrr', (100) 
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where ojk = |fc|. 

We consider an interaction between the electron held and the detector of the form 

= + ( 101 ) 

where x{x),x{x) are currents dehned on the detector Hilbert space JC. 

The stationarity condition, Eq. ([5]), takes the form 


i(x)e-“'|w) = jd\i(g)S^x‘-e‘{q,t))\Lj'), (102) 

x{x)e~'’'^\u) = J- £'(q,t))\oj'), ( 103 ) 


for some \uj') G /C and operators fi^q) and x(9) 

We assume that the detection signal is solely due to the absorption of electrons and not 
due to the absorption of positrons. This implies that the operator xisi) that couples to the 
annihilation operator for positrons annihilates the initial state |a;') of the detector. 

We consider an initial state po for the held that contains N electrons and zero positrons. 
Then, the only correlation functions that contributes into the one-event probability density 
Eq. fl43p are of the form Tr . Hence, all relevant information about the 

initial state is contained in the one-electron reduced density matrix 


pi(A:,r; fc',r') 





( 104 ) 


which is normalized as 


j dfi{k)pi{k,r] k,r) = 1. (105) 

r 

Next, we identify the Hilbert space V of test-functions for the current operators that was 
introduced in Sec. 4.2. V consists of test functions fA on the instantaneous space S, where 
H is a four-spinor index. A general function Ja can be analyzed in terms of the Dirac spinors 
u{k, r) and v{k, r) as 


fA{q) = ^ {e"^''^UA{k,r)4>r{k) -f e^'‘'^VAik,r)x*rik)) , (106) 

where (f>r{k) and Xr{.k) are elements of the Hilbert space 'H 1/2 that carries the irreducible 
representation of the Poincare group for s = |. 

It is important to note that the positive operators T{p) are dehned on the Hilbert space 
V = 7^1/2 © 'Hi /2 rather than on 'H 1/2 as one would expect by a particle treatment of the 
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electron, without taking its QFT coupling into account. Since a spin measurement can only 
have two outcomes, the POVM T(/i) can always be expressed as 


tip) = i (i + ^e) , 


( 107 ) 


where /i = ±1 and S is an operator on V. The specihcation of S and its covariance properties 
depend on the detailed model for the spin measurement, and will be considered in a different 
publication. 

Next, we specialize to the case of a static detector and we express the embedding function 
S(t,Q) = We consider a /i-independent degradation function r], and we evaluate the 

probability density (H3|) . 


PiT,Q,i^) = NC EE d^{k)d^{k')pi{k, r; k', r')uA{k, r)uB{k', r') 

r,r' A,B 

1 


XT] I -{Uk + ujk') ] exp [i{k - k') Q - i{uk - n;fe/)r] 


^AB + I^'Pab^Q-i 2 T A:)) 


(108) 


In deriving Eq. flllOp . we took the limit —)■ 1 and dehned as 'Pab the Weyl-Wigner 

transform of the operator E of Eq. fll07p . 


Sas(Q,P) = j dV"'^(Q-|,i?|T(/x)|g + |,kl) (109) 

We consider the set-up of the ideal time-of-arrival measurement of Sec. 6.2. That is, we 
consider an initial state at the origin of coordinates and a detector localized at x = Q, so 
that we can treat the system as one-dimensional. Treating the position coordinate Q as a 
parameter, and we normalize so that h) coincides with the probability density 

fl66l) . We obtain 


P{t,Q,p) = — EE/ dp,{k)dp,{k')pi{k, r; k', r')uA{k, r)uB{k', r') 


r,r' A,B 


X \/ {Uk + uJk'y - 4m2 exp [i{k - k')Q - i{uk - uJk')T] 


^AB + k^^AB^Q, —{k + k')) 


( 110 ) 


We integrate over r to obtain a probability distribution over p, P{p). Extending the 
range of integration for r to the whole real axis, we obtain 


Pip) = 2 


l + /i^ / dp{k)pi{k,r;k,r')Err'{Q,k) 


where 


P‘rr'{Q,k) = — y^UA{k,r)uB{k,r')T,AB{Q,k). 

2m 


( 111 ) 


( 112 ) 


AB 
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We note that if the Wigner-Weyl transform Tiab is independent of Q, the probabilities 
fillip are also independent of Q. However, this property is not always guaranteed: the 
determination of spin depends on a split of the generator of spatial rotations into an angular 
momentum and a spin part. Since angular momentum depends explicitly on the position 
coordinate, it is expected that the spin operator will also depend on Q. Our analysis strongly 
suggests that the spin operators are contextual, i.e., they depend on the specihc set-up 
through which spin is measured. Further analysis is needed in order to demonstrate this 
point and it will be a subject of a different publication. 

9 Conclusions 

We have elaborated the formalism for relativistic quantum measurements developed in Ref. 
[T] . We constructed a general model for a detector using the master equation approximation 
for the pointer variables, and we identified the mathematical objects necessary for the char¬ 
acterization of a detector. We obtained a simplified expression for the probability associated 
to n measurement events, and we identified regimes that correspond to ideal measurements. 

We applied our results to three classes of systems. In particle detectors, we constructed an 
ideal probability distribution for time-of-arrival measurements and showed that it persists for 
different types of held-apparatus coupling. Then, we constructed the probabilities associated 
to photodetection events. We showed that Glauber’s photodetection theory corresponds to 
the limit of an incoherent detector. Finally, we considered the case of relativistic spin mea¬ 
surements, where we identified a general expression for the associated probabilities involving 
an integral over all possible times of detection. 

In the last section of Ref. [1], we mentioned some possible applications of our formalism 
in relation to the quantum foundations. In addition to this, the elaboration of our results 
on spin measurements and the consideration of photon polarization measurements are direct 
extensions of the work presented here. The strength of this formalism lies in its treatment of 
multiple measurements. The examples that were presented in this paper involved only single 
measurement events, but the extension to n events is straightforward. Of immediate interest 
are (i) the consideration of time-of-arrival correlations in multi-partite systems and their 
relation to entanglement, and (ii) the possibility of using the multi-event probabilities in order 
to provide a Poincare invariant dehnition of quantum information concepts, in particular, of 
qubits and of entanglement. 
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